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O ■ We investigate the finite-temperature spectral functions of heavy quarkonia by 

using the soft- wall AdS/QCD model. We discuss the scalar, the pseudo-scalar, the 
vector, and the axial-vector mesons and compare their qualitative features of the 
melting temperature and growing width. We find that the axial-vector meson melts 

> , 

earlier than the vector meson, while there appears only a slight difference between 

ON . 

| the scalar and pseudo-scalar mesons which also melt earlier than the vector meson. 



I. INTRODUCTION 

The strongly-correlated quark-gluon plasma (sQGP), which is hot and dense matter out of 
quarks and gluons created at the Relativistic Heavy- Ion Collider (RHIC) in the Brookhaven 
National Laboratory (BNL), have been attracting a great deal of interest in its intrinsic 

n n n 

non-perturbative properties [1|, 12J, |3|. Although there is no systematic way to study such 
a non-perturbative system at strong coupling, apowerful technique has developed recently 
based on the gauge/string correspondence [4, a, f|- The idea is that one can treat the 
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strong-coupling regime in the gauge field theory on the boundary by solving the weak- 
coupling string theory (or classical one in the large N c and large 't Hooft coupling limit) in 
the bulk Anti de Sitter (AdS) space. 

A well-known example of successful application of the gauge/string duality to sQGP 
physics is the exact computation of the shear viscosity to the entropy density ratio, i.e. 
rj/s = h/ (ATrks) in an J\f = 4 supersymmetric Yang-Mills plasma 7, s[ 0, [h]]. This value of 
7]/s is much smaller than any observation in reality except for the heavy- ion collisions; the 
hydrodynamic model studies imply that rj/s of QCD matter is as small as suggested by the 



string theory H 



l^ . Besides, it is conjectured that t]/s = h/(4irkB) might be a universal 
lower bound and applied to strong-coupling QCD as well as supersymmetric models. 

The smallness of r]/s is an important indication of the sQGP because a larger reaction 
cross-section leads to a smaller t] in gaseous states. Actually perturbative QCD calculations 
cannot give an account for small r]/s in the weak-coupling regime 13|, [l4j, [l5| . The Monte- 
Carlo simulation of QCD on the lattice is a powerful instrument to look into the non- 
perturbative strong-coupling regime. It is still difficult to estimate t] in fully dynamical 
simulations with quarks, but the (quenched) results so far are not inconsistent with the 
string theory estimate [16|, 1171 ] . 

Another important indication to the sQGP is the in-medium property of heavy quarkonia 
such as J/ip. In the recent lattice QCD simulations the J/ip spectral functions (SPFs) both 
above and below T r have been successfully constructed by means of the maximum entropy 
method (MEM) 18(, which has revealed that the mesonic correlation (a peak in the SPF) 



survives even above twice of T c 19|, |20|, |2l|, |22j, [23|]. It is, however, a non-trivial question 
how to explain such a high melting temperature for J/ip in a conventional way using the 
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28|. We have 



non-relativistic model with the Debye screened potential 
not yet reached a full consensus on the interpretation of the J/ip SPFs above T c , though 
there are many theoretical efforts. Our present aim is to investigate this question using the 
gauge /string duality along the same line as our previous work 29]. 

In Ref. 29[] we calculated the SPFs in the vector channel assuming that the heavy- 
quark sector is decoupled from others. In this work we will extend our analysis to other 
channels; scalar, pseudo-scalar, and axial-vector mesons, namely, Xco, Vc, and \ r -\. Since the 



interpretation of Xco (i- e - whether it melts or not above T c ) was controversial 
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22 



23|, 



it is important to clarify whether cc states in all these channels melt at T > T c and, if so, 
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when they melt. Under the situation that the MEM construction of the SPFs in lattice QCD 
simulations are still difficult for all these channels, it is valuable to take advantage of the 
holographic QCD model to see what spectral shape would transpire in the strongly-coupling 
system. In this work we will use the soft- wall AdS / QCD model 30, |31[ (see Refs. 32|, |33| for 
related works). Although the SPFs at finite temperature and density have been discussed 
by means of the D3/D7 setup 3^], it is not straightforward to deal with the heavy-meson 
SPFs in the D3/D7 model. This is because the only energy scale in this model is fixed by the 
pion decay constant and the SPFs are given as a function of not T and M q independently 
but only T/M q where M q is the quark mass 35|, |36|. The soft- wall model, in contrast, has 
one more phenomenological parameter, c, which is fixed by the meson spectrum. 

One of the important features in the soft-wall AdS/QCD model is that a parameter in 
the bi-fundamental scalar sector controls the chiral symmetry breaking, that is, the chiral 
condensate. Interestingly enough, the chiral condensate is uniquely determined from the 
infrared (IR) boundary condition that is specified by a smooth function in the soft-wall 
model. Therefore, in this model, there is no ambiguity in the IR limit in the case at finite T 



fsee Refs. 
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for holographic approaches to finite-T mesons). The regular IR boundary 
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4l| from which we can 



condition allows us to evaluate the Minkowskian correlator 
compute the SPFs. We will find, in view of the resultant SPFs, that the axial-vector states 
melt faster than the vector ones. We will then clarify its origin in the chiral symmetry 
breaking due to the scalar field which yields a difference between the vector and axial- 
vector equations of motion. On the other hand, the scalar and pseudo-scalar dissociation 
temperatures are almost the same; T ~ T c . Finally, before closing this paper, we shall take 
a closer look at the vector channel. 



II. SOFT- WALL MODEL 

The principle to construct the AdS/QCD model is the bulk/boundary correspondence 
or UV/IR relation. The generating functional in the gauge field theory is equivalent to the 
exponential of an on-shell action in the gravity theory (GKP-W relation 

Z[<f) ] = / e i f dx <t>0(x)O(x)\ _ giSgravityW^ 

\ ' gauge 
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BB), 
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where Z[(f> ] is the generating functional with the source O coupled with an operator 0(x) 
and Cavity is an on-shell action with the boundary condition <p — > <f) at the UV boundary 
where the gauge theory resides. 

The AdS/QCD models are five-dimensional field theories designed to describe QCD prop- 
erties through the bulk/boundary correspondence 32j. The essential ingredients of the 
AdS/QCD model are the AdS space with an IR cutoff (i.e. wall) that translates into a typ- 
ical energy scale in QCD, the Ul(2) x Ur(2) vector fields, Al, A r , and the bi-fundamental 
scalar field X. The vacuum expectation value of X is responsible for the explicit and sponta- 



neous chiral symmetry breaking. The soft-wall model is defined by the following action 3l| ; 



S= / d 5 xe~ cz V=g£ (2) 



£ = tr 



\DX\ 2 + ^ 2 \X\ 2 - -^{g MN g PQ F LMP F L , NQ + g MN g PQ F R , MP F R , NQ ) 



(3) 



where X = X a t a with t a being the generator of £7(2) and D M X = d M X+iA L ^ M X— iXA R ^ M . 
Here M = Xq, Xx, X2, £3, z and g$ are the indices for the five- dimensional coordinates and 



gauge coupling, respectively. We note that g\ = 2An 2 L/N c is concluded by matching |31 ] . 

We also use the Greek index /x = xq,x\,X2,xz to refer to the four-dimensional coordinates. 

The model parameter, c, characterizes the wall location; since the contribution from the IR 

region z > 1/ y/c is suppressed by e~ cz2 , it represents a potential with the wall providing a 

typical QCD scale. The background geometry is specified as the AdS metric as 

7-2 

g MN dx M dx N = — (-dt 2 + dx 2 + dz 2 ) . (4) 

z 

It should be mentioned that the vector meson mass spectra at T = are quantized by 



the normalizability condition and given as the following Regge trajectory 



3l| 



■in 



2 n = 4cn, (5) 



where n is the radial excitation number. Then we can determine c by fitting the above 
relation to the vector meson spectra; p(770), p(1450), p(1700), etc. If we take m p = 0.77 GeV 
for n = 1, we have c = 0.77 2 /4 = 0.148 GeV 2 , while we will later find that the spectral peak 
is slightly shifted from Eq. (jSJ) and will fix c = 0.151 GeV 2 to fit the peak position with the 
mass. 

Now that we fix the model parameter c, let us consider the model at finite tempera- 
ture. Here we shall introduce the following background, which is called the AdS blackhole 
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(AdSBH), 

g MN dx M dx N = (-f(z)dt 2 + dx 2 + —j—dz 2 } , (6) 

with f(z) = 1 — / ' z\ where the horizon is related to the Hawking temperature that is 
interpreted as the system temperature of dual QCD as Zh = l/(nT). It is known that the 
AdSBH is unstable at low temperature, and thus the Hawking- Page-type transition occurs 
at a critical temperature, T c = 0.492a/c. This is a first-order phase transition from the 



42 



43| and can be identified as the color confinement- 



AdSBH to the thermal AdS metric 
deconfinement phase transition. 

In the soft-wall model one can introduce the chiral symmetry breaking explicitly (i.e. 
quark mass) and spontaneously (i.e. chiral condensate) through the bi-fundamental scalar 
field X, which is decomposed as X(x, z) = e 2tn ^ x ' zS) [Xo(z) + S(x, z)\ where x refers to four- 
dimensional coordinates only and Xq is a constant background with respect to x. The 
fluctuations, S and n, represent the scalar and pseudo-scalar fields. In the case when quark 
masses are degenerated, X (z) is proportional to unity in flavor space and satisfies the 
following equation of motion, 

X'^z) + (-2cz + ^)*o(*) + ^7 X °W = °> ( 7 ) 

where the prime stands for the derivative with respect to z. From this differential equation 
we find that in the vicinity of z = the general solution behaves as 

L^ 2 X Q {z)^ l -{M q z+Hz 3 ), (8) 

where, according to the dictionary of bulk/boundary correspondence, the parameters M q 
and S are identified with the quark mass matrix and the chiral condensate, respectively. 
Here, in the soft-wall model, £ is uniquely determined for a given M q so that Eq. ([7]) can 
yield a finite and regular solution of X . This property is a flaw in the light-quark sector 
because M q = always leads to X = and so the spontaneous breaking of chiral symmetry is 



not correctly described unless Eq. ([7]) is modified with higher-order potential terms 44, l45|. 
In the present work, as we discuss later, only the heavy-quark sector is of our interest and 
we need not alter Eq. (JTj) because chiral symmetry is largely broken in an explicit manner. 

We must point out that the conventional soft-wall model has another flaw in the chiral 
properties. In the vicinity of the UV limit two independent solutions of Eq. (JTj) are definitely 
z and z 3 , but if we carefully go beyond the leading order, the former solution receives a 
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correction by a logarithmic term as z - > z + (— 1/2 + log z)z 3 . This higher-order correction 
is small as compared to z, but not small at all to another solution z 3 . Therefore, Eq. ([7]) 
leads to a UV divergent chiral condensate, which is an artifact of the soft-wall model. (There 
is no such logarithmic term in the hard- wall model.) Therefore we need to modify the model 



as done in Refs. 



44 



for example. In the present work we will take the following strategy. 
That is, to solve Eq. ((Tj) numerically, we will force the initial condition by Eq. (jSJ) and find 
S in such a way that the solution in the IR region contains no singularity. This is not a 
fully satisfactory resolution but is acceptable pragmatically for the soft-wall model that is 
only a phenomenological model. 

Let us mention on the asymptotic solutions of Xq near the horizon (z ~ Zh) for the 



finite-T case. Equation (jjj) can simplify by the variable change from z to t = a/3(1 — zJzh), 
which reduces to the Bessel equation near the horizon t ~ 0. It is thus obvious that the 
asymptotic solutions of X are given by the first-kind Bessel function Jo(t) which is regular 
and the second-kind Bessel function Yo(t) which is divergent at t — 0. Since the physical 
solution must yield a finite action, we should pick only Jo(t) up near z — Zh- To this end we 
need adjust an appropriate ratio of M q and E in the initial condition (jSJ) in the UV boundary, 
so that it evolves to Jo(t) near the horizon. We will concretely carry this procedure out in 
later discussions. 



III. FLAVOR-DEPENDENT SOFT- WALL MODEL 

The mass spectra (jSJ) in the soft-wall model successfully reproduce the Regge trajectory of 
the light vector mesons consisting of u and d quarks as seen in the previous section. In order 
to apply this model description to the heavy-quark sector, we propose a modification on the 
soft-wall model in such a way that we treat c as a flavor- dependent parameter. The following 
action in our treatment is composed from two sectors; one is the light-quark (u, d, s) sector 
and the other is the heavy-quark (c) sector; 




where £ heavy takes an almost identical structure with Eq. ([3]) in the light-quark sector. The 
only difference is that fields in £n g ht and in £hea vy belong to U(3) and U(l) groups, respec- 
tively. From the vector-meson mass formula (jSJ), we can determine the model parameters 
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as 

c p = 0.151 GeV 2 , cj/* = 2.43 GeV 2 , (10) 

to reproduce m p = 0.77 GeV and mj/^ = 3.1 GeV (as we have noted, the spectral peak is 
slightly different from Eq. (jHJ) and c is shifted from the naive estimates accordingly). If we 
believe in the mass formula, the above value of cj/^p predicts the mass of the first excited 
state as 4.4 GeV, which overestimates the mass of ip{2S) that is 3.7 GeV. Therefore this 



deviation by around 20% should be taken for a systematic error in this model [38[ . This 
error is of acceptable order as compared to errors associated with other assumptions such 
as the large N c limit, the probe approximation, etc. 

Because of c p <C cj/^, the latter term in the action (jHD is negligible to evaluate the 
magnitude of S. Hence, the critical temperature T c of the Hawking-Page transition in this 
model is solely determined by the former term involving c p , that means T c = 0.492^/c^ = 
0.191 GeV is unchanged. In short, the bulk thermodynamics is dominated by the former 
term, while the heavy-flavor sector is described by the equation of motion deduced from the 
latter term the action ([9]). 

Before closing this section let us comment on possible justification of this model treatment 
with two scales. One may wonder that c should be common to all flavors because it is a 
parameter related to the QCD string tension. Besides, it should be more natural that mm 
arises mostly from M q rather than cj/^p. In the soft-wall model, however, the vector-meson 
field has no direct coupling with Xq and so it does not depend on M q . The important point 
is that cj/^ as a "renormalized" scale can originate from the back- reaction with the heavy 



44 1 



charm-quark mass beyond the probe approximation. In fact it is pointed out in Ref. 
that the back-reaction from the X field in the hard-wall model produces an effective soft- 
wall with c depending on M q . So far there is no such analysis on the back-reaction within 
the framework of the soft-wall model, but it would be a reasonable anticipation that c 
must get larger with heavier M q once the back- reaction is taken into account. The back- 
reaction analysis in the top-down approaches Q, also implies that our treatment could 
be pragmatically acceptable. 



s 



IV. SPECTRAL FUNCTIONS 



In this section we proceed to actual calculation of the SPFs. As seen from the 
bulk/boundary correspondence, we can derive the SPFs in the channel of our interest by 
solving the classical equation of motion in five dimensions. Here we take the AdS radius as 
L = 1 since this quantity disappears in the physical correlation functions. In addition, for 
convenience, we use the dimensionless energy u, momentum q, and temperature t in the unit 
of y/c and change the variable for the fifth coordinate by £ = \/cz, so that we can totally 
eliminate c from the equation of motion. Because only c is a dimensional parameter in the 
model, we can easily restore c to discuss physical quantities. It should be noted that spatial 
and temporal components lead to distinct differential equations since Lorentz symmetry is 
broken in the presence of a medium. For the moment we will focus on the solution of spatial 
fields in this section, then in Sec. IV CI we will address a physics insight into the dependence 
on the polarization direction. 



A. Vector Mesons 



Let us first consider the vector meson whose dual field is Vm = (Al,m + A RjM )/2 and 
then the axial-vector meson whose dual field is Am = (Al^m — Ar,m)/2. We fix the gauge by 
choosing Al, z = A^ z = 0. Besides, we impose 8^Al,h = d^A^ = to get rid of unphysical 
polarization. The linearized equation of motion for the spatial component V x (either x = x±, 
X2, or Xs) of the vector field takes the following form, 



-gg xx g zz (d z v x ) 



+ 



-gg xx d^V x 



0. 



(11) 



Now we move to momentum space by performing the Fourier transformation, V x (x,£) = 
J d A x e Zy ^ p ' x V(p)v(^; p) and substitute the AdSBH metric (\6§ into Eq. (TTTj) . so that we 
reach, 

'3/-. 



v" + 



0. 



(12) 



with = (ou, q 1 , q 2 , q 3 ) and q 2 = {q 1 ) 2 + {q 2 ) 2 + (<Z 3 ) 2 - Here, as we have mentioned before, 
all variables are dimensionless and the prime (') stands for the derivative with respect to £. 

Before solving Eq. ffl2l it would be instructive to pursue the analogy to the Schrodinger 
equation in quantum mechanics. The change of the field, u = (e~ cz2 



-gg xx g zz y/ 2 v, simpli- 
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FIG. 1: Potential U v (£) for the vector fields at dimensionless temperatures, i = T '/ '^/c = 0.05, 0.10, 
and 0.15 at cj 2 = 4 and g 2 = 0. 

fies the equation of motion in the following form; u" — U v (£)u = with the potential, 

^o=« , ^-K 2e+ f)-^ + 0-K7-4 (13) 

Figure [T] shows this potential for various dimensionless temperatures in the unit of yfc. In 
the case at T = (and thus / = 1) the downward-convex potential, £ 2 + 3/(4£ 2 ), yields 
the discrete spectrum, m? = uj 2 — q 2 = 4n (n = 1, 2, . . . ), only for which the wave-function 
is normalizable. We see that the higher t or smaller £/, = l/(7r£) makes the potential less 
convex and eventually it becomes monotonic at t ~ 0.15. With a monotonic potential we 
cannot expect a remnant of the original spectrum any longer. In other words we should 
anticipate dissociation then. Let us confirm this in what follows. 

At finite temperature the potential is no longer rising in the large z side and the normal- 
izability does not quantize the spectrum. We can easily extract the asymptotic solutions of 
Eq. (|12p near the horizon as 

U(0 — C+0+ + C-0- With ± = (l_^)±-&/4 ) (14) 

in the vicinity of £ — > £h- Here <p + represents the out-coming solution and 0_ the in-falling 
solution into the black hole. Near the origin, on the other hand, the solution has the following 
asymptotic form, 

v(0=A$ 1 + B$ , (15) 
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where $1 and $o are two solutions of Eq. ([TT]) satisfying the following UV boundary condi- 
tions; 



7T 

"2 V 



:£Ji(V^m (16) 



a/cu 2 — g 2 

around £ — > 0. Here J\ and Yi are the first-kind and second-kind Bessel functions, re- 
spectively. In the above we normalized <3>i and <3>o m such a way that $i(£ = e) = 1 and 
$o(£ = e) = e 2 and also we assumed that uj 2 > q 2 . In the case that q 2 > uo 2 we should replace 
the above by $ x -> v^-cu 2 ^^ 2 -^) and $ -> (2/y/q 2 -u 2 ) £ h(^/uj 2 -q 2 £)• In 
what follows we will fix the overall normalization of v(£) by adopting the commonly used 
prescription, A = 1, so that B should be unique once the IR boundary condition is specified. 
Following the procedure elucidated in great details in Refs. [39|, |40j, |41] we can compute 



the Green's function in Minkowskian space-time. The IR boundary condition must be 
v(£ — > £h) — c„0_ (i.e. c + = 0) to acquire the retarded Green's function according to 



Ref. [39j. We can make v(£) satisfy this IR boundary condition by choosing B appropriately 
at £ ~ (where A = 1 is chosen so that v(£ — > e) = 1); Then, B, which is now a complex 
number, is uniquely fixed by the IR boundary condition; 



v(0 = ^i(0 + B(u,q)%(0 



as 



(17) 



As we defined above, we can generally solve $i and $o from Eq. (TTTj) from the UV asymptotic 
forms fTT6l) toward the IR side. If we have, 



*<(£) 



(18) 



where i = 0, 1, then we can readily conclude B(lj, q) = —ai/a . 

Once B is obtained, the bulk/boundary correspondence ([TJ allows us to compute the 
Green's function, that is given as 



D R (cu, q) = -CBmQuV 



-2C 



- ln^v^^e) 



(19) 



where C is a constant given as C = N 2 / (6Air 2 L) . The spectral function is, by its definition, 

(20) 

Here we note that only B(u,q) has an imaginary part in Eq. (fT9~j) . 



1 2C 
p(uj,q) = lmD R (uj, q) = — lmB(uj,q). 

IX 71 
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FIG. 2: Vector spectral functions ImB(u>,0) for the temperatures, t = 0.07, 0.09, 0.10, and 0.13. 
The melting temperature for the lowest-lying peak is about t = 0.14. 

We are now ready to come to the numerical calculation. We plot lmB(u, q) by calculating 
a,i(u), q) numerically as a function of u and q at various temperatures and show the SPFs at 
q = in Fig. [2J We should remark our convention that we refer to lmB(u>,q) as the SPF 
neglecting an overall factor. 

B. Axial- vector Mesons 

Next we shall move to the SPFs in the axial-vector channel. We can follow exactly the 
same procedure as the previous one to look into the axial-vector fields, which we denote as 
Am = (Al,m — Ar,m)/2- We again define the Fourier mode of the transverse component 
A x (where either x = Xi, £2, or 23), i.e. a(£;p). The dimensionless equation of motion is 
expressed as 

a"(0 + (^p - 2*) a'(0 + " 7) «(0 + ^7 X o «*) = 0- (21) 

We see that the above (1211) is just the same as Eq. (1111) for the vector fields except for the 
last term involving Xq where X® is a solution of Eq. (CO). The chiral symmetry breaking from 
M q ^ and E ^ is introduced by X ~ \{M q £ + ££ 3 ) near £ = and is responsible for 
the mass splitting between the vector and axial-vector channels. We can also expect the last 
term becomes negligible as compared to the third term for large u or q, so that the highly 
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excited radial states exhibit degeneracy between the vector and axial- vector mesons 45], 



481 ] . As discussed in Sec. flU the 



which has been observed in the excited baryon spectrum 
quark mass M q and the chiral condensate £ are not independent in the soft-wall model and 
once M q is fixed, £ is uniquely determined so as to yield a regular solution of Xq in the IR 
region under a requirement that the UV initial condition is forced to be Eq. (jSJ). 
Now we fix M q as the charm mass; 

M q = M charm = 1.30 GeV, (22) 

to derive the associated chiral condensate £ by the shooting method numerically. In our 
calculation we obtain £ ~ —(3.1 GeV) 3 , which seems overestimation but within a reasonable 
range of order. Using these M q and £ we can get a regular numerical solution of the 
background scalar field X (^). One noticeable fact to be mentioned is that, since X Q is 
regular for an appropriate choice of M q and £ both near the horizon £ ~ £/, and near the 
boundary £ = 0, the boundary conditions for a(£) are (almost) the same as those for v (£) 
as follows; 

a(£) — ► c + + + c„0_, (23) 
near £ — > and near the UV boundary we have 

a(0 = + B$' 0J (24) 

where $' x and $' are two solutions of Eq. ( I2T1) satisfying Eq. ( fl6l) with cj 2 — q 2 replaced by 
oo 2 - q 2 + (24vr 2 /N C )M 2 . 

Hereafter we will trace the same analysis from Eq. ([TBI to Eq. (1201) in the previous 
subsection. In the picture of the Schrodinger equation the corresponding potential for the 
axial vector case is given by 

U (!) - S 2 + 3 r (/ ' )2 + f " 1 a 2 ) + 96712 X 2 (25) 

Ua^-Z + ^-J^ + l)-^p + 2f-j{j- q ) + N^f X ^ (25) 

We show the profile of U a (£) in Fig. (3) It is clear in view of Figs. [1] and [3] that the axial- 
vector potential becomes less downward-convex earlier than the vector case, and the shape 
looks monotonic already around t ~ 0.10. Thus we can anticipate that the axial-vector 
spectral peaks should melt much earlier than the vector ones. In fact t ~ 0.10 corresponds 
to T = 0.10-y/c = 0.16 GeV, which is below the deconfinement temperature T c = 0.191 GeV, 
meaning that the axial- vector mesons should melt at the phase transition. 
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FIG. 3: Potential U a (£) for the axial-vector fields at dimensionless temperatures, t = 0.062, 0.07, 
and 0.10 at u? = 4 and q 2 = 0. 

Now let us derive the axial-vector SPFs. The solution satisfying the in-falling boundary 
condition determines a complex value of B(uj,q), 

a(£) = & 1 {Z) + B(u,q)& (t) — c_<M0 as £ - (26) 

Through the same procedure as elaborated in the previous subsection, we estimate the SPFs 
for the axial vector mesons by evaluating lmB(u, q) numerically. We show our numerical 
results for the axial-vector SPFs with q = in Fig. HI Here we depict lmB(u, q) divided by 
10 to make its scale similar to Fig. [2j The overall factor takes a different value depending 
on the vector and axial- vector channels because of our normalization convention A — 1. 
Therefore, under the choice of A = 1, it is not a physically meaningful comparison to take 
the absolute magnitude of spectral heights seriously. 

C. Discussions for Vector and Axial- vector Mesons 

Here let us make a comparison between the vector (Fig. [2]) and axial-vector (Fig. Hj) 
channels. For low temperatures, the lowest-lying peaks are located at u 2 = 3.92 for the 
vector case and to 2 = 4.72 for the axial- vector case. If we fix c = 2.43 GeV 2 to reproduce 
m J/ip = 3.1 GeV, then we have the mass in the axial- vector channel as m Xcl = 3.4 GeV, 
that is in agreement with the experimental value 3.51 GeV. 
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FIG. 4: Axial- vector spectral functions Im5(w,0) x 10 _1 for the temperatures, t = 0.062, 0.07, 
0.08, and 0.10. The melting temperature for the lowest-lying peak is about t ~ 0.10. 

It is an apparent feature seen in Figs. [2] and H] that the spectral peaks become more 
collapsed and the peak positions move smaller as t increases. We also note that the second 
lowest-lying states melt far earlier than the lowest-lying states both in the vector and axial- 
vector cases. This is quite natural because higher excited states are less stable generally. In 
terms of the potentials illustrated in Figs. [Hand [3] a larger uo 2 causes stronger absorption into 
the black hole by the term, —u 2 /f 2 , which is negative large near the horizon. Furthermore, 
as seen from the t = 0.07 curve in Fig. [2} the lowest-lying state moves smaller only slightly, 
while the excited states shift more drastically. These qualitative properties of the SPFs are 
consistent with the lattice QCD observations for the heavy quarkonia. 

As seen from Figs. [2]and|4l the dissociation takes place around T ~ 0.15a/c ~ 0.23 GeV 
for the vector lowest-lying peak and around T ~ 0.10-^/5 ~ 0.16 GeV for the axial- vector 
one, where c = cj/^ = 2.43 GeV 2 as discussed before. The deconfmement transition occurs 
at T c = 0.492a/c = 0.19 GeV as mentioned in Sec. IHIl Thus, in our soft-wall QCD model, 
the vector charmonium J/ip melts above the critical temperature; T ~ 1.2T C , while the 
dissociation temperature of the axial-vector charmonium Xci is much lower; T ~ 0.8T C , 
which indicates that xa does not survive above T c and melts suddenly at the deconfmement 
transition. It is obvious in our argument that the chiral symmetry breaking induced by Xq 
causes this difference between the vector and axial-vector SPFs. 
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D. Scalar and Pseudo-scalar Mesons 

Here we go on to the SPF for the scalar and pseudo-scalar fields, whose lowest-lying peak 
can be identified as Xco and r] C Q. As discussed in Sec. El we can introduce the dual fields 
of the scalar and pseudo-scalar mesons by decomposing the bi-fundamental scalar field as 
X = e 2lU ^ x ' z ' > [X (z) + S(x, z)], where X is the background part, S is the scalar field and IT 
the pseudo-scalar field. We will denote the Fourier modes of S and II as s and tt, respectively, 
in what follows below. 

Before addressing the SPFs, we need to consider the holographic renormalization and 



counter terms 



,|50|, Q to give regular results near the boundary as well as physically 
meaningful SPFs in the scalar channel. The action with respect to the scalar and pseudo- 
scalar fields in the quadratic order of S 2 and IT 2 is given by 

S= / d A x^-(-S'S-4X^U'U) +S com (27) 

Jz=0 Z 

after the integration by parts, in which the functional derivative of S eom leads to the equations 
of motion. The first term is UV divergent at z — > and requires the renormalization counter 
term that is constructed in such a way that the covariance holds; 

Sr C n= I d 4 xe~ cz2 ^X 2 , (28) 



2 = 



2 ~,-2\ 



Z ,Z ,Z 



49 



50 



5l|. 



where 7 is the determinant of the induced metric defined as 7^ = diag(— fz~ 2 
and thus -^—7 ~ z~ 4 near z = 0. The renormalized action is defined as S + Sj 
Once we comply with this renormalized procedure, we can follow the same procedure as in the 
previous case for the vector and axial-vector mesons. In this way we find the dimensionless 
equation of motion for the scalar and pseudo-scalar fields; 

s » K) + (/_l_ 2f ), (f) + (^_^ + _L) sK) = , (29) 

-"(«) + (^ - 2f + ^) AO + (j 2 - y)"({) = 0. (30) 

Here the dependence on the quark mass and the chiral condensate is introduced into the 
pseudo-scalar solution through the second (first-derivative) term in Eq. ( 1301) . while the scalar 
equation of motion does not have such a term. This difference should be attributed to distinct 
mass spectra and dissociation temperatures between the scalar and pseudo-scalar mesons 
like the vector and axial-vector cases. 
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FIG. 5: Potential U s (£,) for the scalar fields at dimensionless temperatures, i = 0.05, 0.10, 0.12, 
and 0.15 at J 1 = 4 and g 2 = 0. 
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FIG. 6: Potential U n (£) for the pseudo-scalar fields at dimensionless temperatures, t = 0.05, 0.10, 
and 0.15 at J 1 = 4 and g 2 = 0. 
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Here again, we shall trace the same procedures as those from Eq. ffT3l) to Eq. (|20|) . The 
potentials for the scalar and pseudo-scalar fields in the picture of the Schrodinger equation 
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are read from the equations of motion as 



^2 , 15 /' , 3 



Us(o = e + ^- 2 -^(2z+-^)-^ + ^-(^--), 



rr ,^ ,2 15 /' / , 3\ (/') 2 /" /u; 2 fc 2 \ 3 



(f) 2 


+ r_ 

2/ 




A; 2 


4/2 


U 2 


" / 


(f) 2 


+ £- 
2/ 




k 2 


4/ 2 


^/ 2 


" / 



4£ 2 2/ v ey 4/ 2 2/ v/ 2 /ye 2 / 

Remarkably, the potential for the pseudo-scalar field is independent of the background solu- 
tion X since all the terms depending on X are put together into a form of the equation of 
motion (J7|). The difference between the scalar and pseudo-scalar fields is only the last term 
in Eq. (1321) . We depict these potentials in Figs. OandO The results look very similar and 
turn monotonic around t ~ 0.12. 

Then we find that the behavior of the solutions in the near-horizon region is again given 
by 0±- I n the opposite side of the UV limit we have two solutions for the scalar and pseudo- 
scalar fields. That is, 

s(0 = A a &[ + B S §1 tt(0 = + B^ , (33) 

where $o an d $i are defined in Eq. f|T6|) and $q and are the solutions of the equation of 
motion ( |29l) with the boundary conditions; 3>q(£ = e ) = 6 an d = e ) = e 3 - 

Here let us note that, strictly speaking, the scalar field corresponding to the scalar source 
at the boundary is s(£)/£, and thus the boundary solutions behave asymptotically as 1 and 
£ 2 like the other channels. Therefore the SPFs are characterized by the imaginary part of 
the complex coefficients B s and B w . To calculate the retarded Green's function we fix B s 
and B n requiring the in-falling boundary condition near the horizon. 

We numerically calculate lmB s and Imi^ and make plots as a function of u at q = 
in Figs. [7] and El Here we present Imi?,,- divided by 10 3 since the normalization A = 1 
gives an irrelevant overall factor again. The lowest-lying peak at small temperature for the 
scalar and pseudo-scalar channels is located at uj 2 = 5.85. The model outputs for the meson 
masses are thus m Xc0 m m Vc = 3.8 GeV, which is not good as compared with m Vc = 3.0 GeV 
and m Xc0 = 3.4 GeV. The scalar spectral peaks are located around uj 2 = {An + 6)c with 



52]. Here we 



n — 0, 1, 2, . . . for low temperatures, as is consistent with the results in Ref. 
see that the pseudo-scalar spectral peaks are found at nearly the same positions as the scalar 
mesons. 

The lowest-lying peaks are gradually collapsed and moving smaller as the tempera- 
ture increases, while the excited peaks dissociate much earlier and shift more drasti- 
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cally. These spectral patterns are qualitatively similar to the vector and axial-vector 
cases. We can observe that the lowest-lying spectral peaks melt out around t ~ 0.13, 

1. e. T = 0.13a/c = 0.20 GeV for the scalar and pseudo-scalar channels both, which is slightly 
above the deconfinement temperature; T ~ 1.05T C . If we take a closer look at the respective 
SPFs, we notice that the scalar meson melts only slightly earlier than the pseudo-scalar 
meson. The difference is, however, hardly perceivable and we can say that the scalar and 
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FIG. 9: Mass shift squared as a function of the width with changing temperatures. The dashed 
curve smoothly connects the calculated points at various temperatures for the vector and axial- 
vector channels. 

pseudo-scalar channels are degenerate regardless of the chiral symmetry breaking. 

V. MORE DISCUSSIONS ON THE VECTOR MESONS 

We have seen that only the vector meson, i.e. J/ip, survives above T c (up to T ~ 1.2T C 
in our model), the axial- vector Xci suddenly disappears at T = T c , and the scalar Xco an d 
pseudo-scalar rj c immediately melt around T ~ 1.05T C . Therefore, it should be worth while 
taking a more serious look at the vector SPFs only. In this section we analyze the vector 
SPFs by deducing the relation between the mass shift Am and the width broadening T with 
changing t. We also discuss the evolution of the SPFs at finite momentum q. Finally we 
briefly mention on the dependence on the polarization direction. 

A. Mass Shift and Width Broadening 
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According to our previous work 29], a functional ansatz, au h /\{uj — ujq) 2 + T 2 ] can fit 
the SPFs pretty well. In this way we can numerically read the peak position Uo(t) (leading 
to the mass shift defined by Am(t) = ujq(0) — U)o(t)) and the width T(t) determined as a 
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function of t. 

In Fig. [9] we plot the mass shift squared (Am) 2 and the width T associated with the 
lowest-lying peak in the vector and axial- vector SPFs. It is an intriguing finding from Fig. [9] 
that, even though the SPFs shown in Figs. [2] and H] look similar at a glance, the qualitative 
behavior of the mass shift is completely different. In the case of the vector meson Am 2 is 
saturated as t increases, while T continues growing. The relation between Am and T has 



been investigated in the QCD sum rule 



531 ]. which is seemingly inconsistent with Fig. [9] 



but a careful consideration clarifies consistency [54J. In the previous work in Ref. [29j we 
proposed a definition for the dissociation temperature by means of the saturating behavior 
of Am 2 around t = 0.14. This working definition works for the vector meson, while the 
axial-vector peak keeps becoming lighter (i.e. larger Am 2 ) and broader (i.e. larger T) and 
thus there is no saturation observed. It is an interesting question whether our prediction 
about the relation between Am 2 and T in the axial-vector channel can be confirmed or not 
in other models such as the QCD sum rule. 

B. Finite Momentum 

In this subsection we briefly discuss the momentum dependence of the SPFs. There are 



several lattice QCD results on the J / ip SPFs at finite momentum |55l . l56| . Although it is 
not clear whether the lattice simulation achieves accuracy enough to be reliable, the general 
tendency is that the spectral peaks are attenuated as q get larger. 

Here we present the results only for the vector channel because only J/ ip survives above 
T c in the soft-wall model, which is our finding in this paper. We plot the numerical results 
in Fig. [10] for q 2 ranging from to 12 with t = 0.10 fixed. We choose this temperature to 
make it easier to grasp the qualitative feature of the fairly prominent peaks in the SPFs, 
though we know that t = 0.10 is below T c . The conclusion is, of course, unaltered even if 
we carry the analysis out on the case at T > T c as long as the peak remains. 

It is apparent that the spectral peak is gradually collapsed as q increases. This result 
is quite non-trivial and peculiar to the non-perturbative regime since in the perturbative 



evaluation a larger q makes the spectral peak less sensitive to the medium effect 571] - It has 



been studied in Refs. 



58 



] that, in the strongly-correlated N = 4 Super Yang-Mills theory, 



J/i/j melts at high q, or in a frame where J/ip is at rest, it melts under the hot wind of 
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FIG. 10: Spectral functions lmB(u),q) as a function of to and q for a fixed temperature t = 0.10, 
where u, q, and t are all dimensionless in the unit of \fc. 



QGP matter. This conclusion has been confirmed in the top-down holographic QCD model 



later 



59] . The discussions in Refs. 58|, |59j did not originate from the shape of the SPFs, 



however. Our present results add another confirmation of the hot screening scenario, and 
maybe the first evidence directly inferred from the shape of the SPFs. 



C. Polarization Dependence 

We finally discuss the vector solution of the equation of motion ffl"2]) . We can easily 
see that the equation of motion takes slightly different forms depending on whether the 
polarization is V x or Vq. One may think that this difference should be interpreted as distinct 
behavior of the transverse and longitudinal modes in a medium. In fact, usually, if we see 
some vector-like collective excitations with finite momentum q that is directed to q$ for the 
moment, a linear combination of 0th and 3rd components describes the longitudinal mode, 
which becomes distinguishable from the transverse modes. 

In this sense, it is surprising that Eq. (fl~2l) takes completely the same form regardless of 
the choice of x — x\, x 2 , or x 3 , that is, Eq. (fT2|) is insensitive to whether is parallel or 
perpendicular to q. This is a very strong statement. Usually, the transverse and longitudinal 
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polarizations become degenerated only when q = because the rotational symmetry is 



restored then 



60j. (See also Ref. 6l| for the hydrodynamic limit of the longitudinal mode). 



We point out that the above-mentioned statement can translate into the interpretation 
that there is no jet in the strong coupling regime 62|. In our case the equation of motion 
and thus all the physical results are given as a function of q 2 only and V^'s are completely 
equivalent for and £3. This means that the rotational symmetry of the system 

is not broken at all even though a finite momentum q is inserted. One can intuitively 
understand this as quick realization of the equipartition of inserted momentum over the 
system. Such a picture is very similar to the finding of Ref. 62] that no jet can exist in 
a strong-coupling M = 4 supersymmetric Yang-Mills medium. Actually, if we assume the 
vector dominance, we can apply our results of the vector SPFs directly for the problem of the 
dilepton production, which may be an interesting direction regarding the future extension 
of our work. 



VI. SUMMARY 

In this paper we derived the SPFs of meson states in the vector, axial-vector, scalar, 
and pseudo-scalar channels at finite temperature using the soft- wall AdS/QCD model. We 
pointed out that the SPFs in these channels have several qualitative features as follows: 

(i) Only J/ip survives above the deconfinement transition up to T ~ 1.2T C and x& 
completely melts at the transition. The scalar Xco and pseudo-scalar rj c are almost degenerate 
in our model and melt soon above T c . 

(ii) The relation between the mass shift squared Am 2 and the width T is characteristic 
to J/ip and Xd- in the vector channel Am 2 looks almost linearly proportional to T at 
small temperatures until it is saturated at the dissociation. In the axial-vector channel, in 
contrast, both Am 2 and T keep growing up with increasing temperature. 

(iii) The spectral peaks diminish at finite momentum, as is consistent with the scenario 
of the J/ip suppression under a hot wind of QGP matter. 

(iv) All the results on the vector and axial- vector meson properties respect the rotational 
symmetry regardless of the presence of the momentum insertion. This should be interpreted 
as the equipartition of the momentum in a medium in the strong-coupling regime. 
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For more realistic studies to investigate the non-perturbative aspect of QCD, we need 
to construct a better model than the soft-wall QCD model that we adopted in this work. 
In the process of concrete computations, in fact, we realized that the conventional soft- 
wall model does not satisfy the requirement that the bi-fundamental scalar field should be 
Xq ~ M q z + £z 3 near the UV boundary (z ~ 0) but leads to a logarithmic correction 
z 3 log(z). The presence of z 3 \og(z) in the solutions of the equation of motion brings huge 
uncertainty in evaluating the chiral condensate £ numerically for a given quark mass M q . 

In addition to this problem of the asymptotic solution, there is another problem, that 
is, the conventional soft-wall model cannot describe the chiral phase transition. In reality 
what should be expected is that chiral symmetry is restored at the deconfinement transition 
simultaneously and then the vector and axial-vector channels become identical. In our case 
the quark mass is significantly heavy and breaks chiral symmetry badly, and thus we can 
consider that the lack of chiral restoration does not affect our results. Nevertheless, it is not 
clear a priori if not only M q but also E have a substantial effect on the discrepancy between 
the vector and axial-vector mesons. To circumvent all these problems we will be able to use 



the modified soft-wall model 
Sakai-Sugimoto models 631 ] . 
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451 ] or the top-down approaches such as the D3/D7 and 



There are many directions in which the present work can be extended in the future. One 
example is the application to the dilepton production problem for which the vector SPF is 
the essential ingredient. We could maybe use more realistic holographic models mimicking 
the QCD equation of state 64|. It is also an interesting generalization to introduce not only 
the temperature effect but also the baryon density or the baryon chemical potential. Then, 
the Chern-Simons coupling mixes the vector and axial- vector mesons 65J, which leads to an 



additional spectral broadening 



66|]. 
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